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Background: R-process nucleosynthesis models rely, by necessity, on nuclear structure models for input. Partic¬ 
ularly important are beta-decay half-lives of neutron rich nuclei. At present only a single systematic calculation 
exists that provides values for all relevant nuclei making it difficult to test the sensitivity of nucleosynthesis mod¬ 
els to this input. Additionally, even though there are indications that their contribution may be significant, the 
impact of first-forbidden transitions on decay rates has not been systematically studied within a consistent model. 

Purpose: To provide a table of /3-decay half-lives and /3-delayed neutron emission probabilities, including first- 
forbidden transitions, calculated within a fully self-consistent microscopic theoretical framework. The results are 
used in an r-process nucleosynthesis calculation to asses the sensitivity of heavy element nucleosynthesis to weak 
interaction reaction rates. 

Method: We use a fully self-consistent covariant density functional theory (CDFT) framework. The ground state 
of all nuclei is calculated with the relativistic Hartree-Bogoliubov (RHB) model, and excited states are obtained 
within the proton-neutron relativistic quasiparticle phase approximation (pn-RQRPA). 

Results: The /3-decay half-lives, /1-delayed neutron emission probabilities and the average number of emitted 
neutrons have been calculated for 5409 nuclei in the neutron-rich region of the nuclear chart. We observe a 
significant contribution of the first-forbidden transitions to the total decay rate in nuclei far from the valley of 
stability. The experimental half-lives are in general well reproduced, both for even-even, odd-A and odd-odd 
nuclei, in particular for short-lived nuclei. 

Conclusions: In certain regions of the nuclear chart, first-forbidden transitions constitute a large fraction of 
the total decay rate and must be taken into account consistently in modern evaluations of half-lives. Both the 
beta-decay half-lives, and beta-delayed neutron emission probabilities have a noticeable impact on the results of 
heavy element nucleosynthesis models. 

PACS numbers: 21.10.Tg, 21.60.Jz, 23.40.-s, 26.30.Hj 


I. INTRODUCTION 

One of the currently most active fields of nuclear astro¬ 
physics deals with the synthesis of elements heavier than 
iron by the r-process. Even though its astrophysical site 
has not yet been identified, it is commonly accepted that 
it occurs in explosive environments involving relatively 
high temperatures (up to 1 GK) and neutron densities 
(n n > 10 2 ° g cm -3 ). Apart from the complexities in¬ 
volved in the astrophysical modeling, the r-process rep¬ 
resents a particularly difficult challenge due to the large 
amount of nuclear input required. It is a complex, dy¬ 
namical process involving a delicate interplay between 
the strong, electromagnetic and weak interactions, re¬ 
quiring the knowledge of a number of observables in sev¬ 
eral thousand nuclei across the whole nuclear chart. One 
of the crucial nuclear properties having a direct impact on 
the distribution of elemental abundances is the /3-decay 
half-lives of the participating nuclei [1, 2]. It determines 
the speed of matter flow towards higher atomic numbers, 
setting the time scale for the r-process. /3-delayed neu¬ 
tron emission occurs during the whole r-process duration 
but becomes particularly important at late phases when 
a competition between neutron captures and beta-decays 
takes place [3]. Particularly important are the beta- 


decay half-lives of nuclei around magic neutron numbers 
N = 50, 82 and 126. The neutron separation energies 
show discontinuities around these magic numbers result¬ 
ing in rather low neutron capture rates. As a consequence 
the r-process matter flow moves closer to the stability 
where nuclei have substantially larger beta-decay half- 
lives. Thus matter accumulates around magic numbers 
N = 50, 82 and 126 producing the observed peaks in the 
solar system r-process distribution. 

At present only a few half-lives of r-process nuclei in 
the vicinity of the magic numbers N = 50 [4-7] and 
N = 82 [8-11] have been measured. Even though the ex¬ 
periments at GSI have provided valuable information for 
half-lives approaching the third r-process peak [12-16], so 
far no experimental half-lives for r-process nuclei at the 
N = 126 shell closure are available. Hence, r-process nu¬ 
cleosynthesis calculations rely mainly on theoretical half- 
lives. The interacting shell model [17] has recently been 
extended to include first-forbidden transitions and has 
been applied to the calculation of half-lives for r-process 
nuclei around N = 50, 82 and 126 shell closures [18]. 
However, due to the increasing computational cost, shell- 
model calculations are restricted to nuclei near closed 
shells. For this reason most of the half-lives calculations 
for r-process nuclei are based on the quasiparticle random 
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phase approximation on top of semiempirical global mod¬ 
els [19-21], the Hartree-Fock-Bogoliubov [22], or the rel¬ 
ativistic Hartree-Bogoliubov [23-25]. The quasiparticle 
random phase approximation (QRPA) approach based on 
a global effective interaction presents a valid alternative 
to the interacting shell model. It has been successfully 
employed in numerous applications [26] and can provide a 
systematic description of /3-decay properties of arbitrar¬ 
ily heavy nuclei. The importance of performing calcula¬ 
tions based on the self-consistent mean-field models, as 
opposed to the empirical mean-field potentials has been 
already emphasized [20]. Currently, however, the only 
available large-scale calculation of /3-decay properties of 
nuclei used in heavy element nucleosynthesis simulations 
is based on an QRPA calculation based on a schematic in¬ 
teraction on top of the Finite Range Droplet Model [19]. 
Additionally, in /3-decay calculations the impact of first- 
forbidden transitions on the total decay rate is mostly 
unexplored. While these transitions and their higher or¬ 
der counterparts have been included in the investigations 
of electron and muon capture rates and neutrino-nucleus 
interactions, their use in /3-decay studies has been either 
limited to particular isotopic chains [27] or have been de¬ 
scribed using a different theoretical foundation such as 
the gross theory [19]. Furthermore, due to the existence 
of a single data table of /3-decay half-lives, it is not pos¬ 
sible to perform meaningful sensitivity studies on the in¬ 
fluence of beta-decay rates on r-process nucleosynthesis. 

In an effort to avoid effects coming from combining 
models built on differing theoretical bases, and with the 
aim of providing a new, high precision data table of /3- 
decay properties, in this study we utilize a fully self- 
consistent theoretical framework based on the relativistic 
Rartree-Bogoliubov (RBB) model [28] for the descrip¬ 
tion of the ground state of open- and closed-shell nuclei 
with the proton-neutron relativistic QRPA (pn-RQRPA) 
where the residual interaction is derived from the same 
density functional as was used for the ground state cal¬ 
culations. In this way, we ensure the consistency among 
the various required properties, from the masses and Qp 
values to the nuclear response strength functions. The 
RBB + pn-RQRPA was already successfully employed 
in the study of Gamow-Teller and higher order reso¬ 
nances [29, 30], /3-decay half-lives [23, 24] and stellar 
weak-interaction processes [31, 32]. This framework also 
enables the treatment of first-forbidden transitions on an 
equal footing with the Gamow-Teller transitions, allow¬ 
ing for a meaningful investigation of the contribution of 
the J > 0 modes of nuclear excitation and for studying 
the regions of interest where this contribution is signif¬ 
icant. We use this model to calculate the /3-decay half- 
lives and the ^-delayed neutron emission probabilities of 
5409 neutron-rich nuclei with 8 < Z < 124, including 
even-nuclei but also nuclei with an odd number of parti¬ 
cles. In Section II we present the theoretical formalism 
of the model, while Section III contains all the neces¬ 
sary expressions for the evaluation of matrix elements 
and nuclear properties, together with the presentation 


and discussion of obtained results. Finally, in Section IV 
we provide the concluding remarks. 


II. THEORETICAL FORMALISM 

A. QRPA calculations 

The calculation of beta-decay rates requires the cal¬ 
culation of both the nuclear ground state and excited 
states of the daughter nucleus and the transitions be¬ 
tween them, together with the evaluation of the lep¬ 
ton phase space involved in the transition. Because 
the calculation requires a good description of a wide 
range of physical quantities, and because the goal is 
to obtain the decay rates for a very large range of nu¬ 
clei, we employ a fully microscopic theoretical frame¬ 
work based on the relativistic nuclear energy density 
functional (RNEDF). The RNEDF based framework em¬ 
ploys the self-consistent mean field for nucleons and a 
minimal set of meson fields; the isoscalar scalar cr me¬ 
son (J 71 ' = 0 + ,T = 0), the isoscalar vector oj meson 
(J 77 = 1 ~,T = 0) and the isovector vector p meson 
(J 77 = 1~ ,T = 1), supplemented with the electromag¬ 
netic field. The meson-nucleon interaction is included 
with a minimal set of the interaction terms, where the 
vertex functionals include explicit dependence on the 
nucleon vector density [28]. The nuclear ground state 
properties are described using the relativistic Rartree- 
Bogoliubov model (RBB), which properly describes the 
pairing effects in open-shell nuclei. For the model param¬ 
eters that determine the density-dependent couplings and 
the meson masses, in this work the D3C* parametrization 
is used [24], as it was previously shown that it provides 
a good description of /?-decay half-lives in medium and 
heavy nuclei. The pairing correlations in open shell nu¬ 
clei are described by the finite range Gogny interaction, 
with the D1S parametrization [33]. 

The excited states are obtained using the proton- 
neutron relativistic quasiparticle random phase approxi¬ 
mation (pn-RQRPA), formulated in the canonical single¬ 
nucleon basis of the RBB model [34] and extended 
to the description of charge-exchange excitations (pn- 
RQRPA) [29]. The RBB + RQRPA model is fully self- 
consistent in both the ph and pp channels. The same 
interactions are used in the RBB equations that deter¬ 
mine the canonical quasiparticle basis, and in the matrix 
equations of the RQRPA. Transitions between the 0 + 
ground state of a spherical parent nucleus and the J 77 
excited state of the corresponding daughter nucleus are 
induced by a charge-exchange operator T JM . Assuming 
spherical symmetry of the nuclear system, the quasipar¬ 
ticle pairs can be coupled to good angular momentum 
and the matrix equations of the pn-RQRPA read [35]: 
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( A J B J \ ( X xj \ _ ( 1 0 \ ( X xj \ 

^ B J A‘ J ) ^ F AJ ) ~ ^ 0 -1 ) ^ Y XJ ) ’ 

(!) 

where the matrices A and B are defined in the canonical 
basis [35]. For each energy E\, X xj and Y XJ in Eq. (1) 
denote the corresponding forward- and backward-going 
QRPA amplitudes, respectively. The transition matrix 
element between the ground state of the parent nucleus 
and the excited state of daughter nucleus, induced by the 
operator T JM , reads 
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B 


± _ 
A J ~ 


'52(p\\T J \\’n) (Xp^UpVn + (-1 ) J Y™v p u n ) 


( 2 ) 

where the X and Y amplitudes are obtained from solving 
the pn-RQPA equation (1). 

In the T = 1 channel of the pn-RQRPA we use the 
pairing part of the Gogny force: 


V PP {1, 2 )=Y1 (Wi + BiP° - HiP T 

2 = 1,2 

—MiP a P T ), (3) 


with r 12 = |»T ^ 7*21, P a ,P T the spin and isospin ex¬ 
change operators and the D1S set of parameters m, Wi, 
Bi , Hi and Mi (i = 1,2) taken from [33] . This force 
has been carefully adjusted to pairing properties of fi¬ 
nite nuclei all over the periodic table. In particular, the 
basic advantage of the Gogny force is the finite range, 
which automatically guarantees a proper cut-off in the 
momentum space. For the T = 0 proton-neutron pairing 
interaction in open shell nuclei we use a form consist¬ 
ing of a short-range repulsive Gaussian combined with a 
weaker longer-range attractive Gaussian. 

2 r 2 

V 12 = — Vo'^gje V fls = i i 7’=o- (4) 

i=i 


where lis=i,T=o projects onto states with 5=1 and 
T = 0. The ranges /.q = 1.2 fm and g 2 = 0.7 fm of 
the two Gaussians are taken from the Gogny interaction. 
The relative strengths g i = 1 and g 2 = — 2 are chosen 
so that the force is repulsive at small distances. The 
only remaining free parameter is Vo, the overall strength. 
This interaction, with a constant value of Vo was used 
in the non-relativistic QRPA calculation [22] of /3-decay 
rates for spherical neutron-rich r-process waiting-point 
nuclei. Two relativistic calculations of /3-decay half-lives 
of neutron-rich nuclei [23, 24] have shown that a single 
value of the overall interaction strength cannot be suc¬ 
cessfully used in different regions of the nuclear chart. 
Thus, we take the ansatz proposed in [25] 

Vo = Vl+ ^ _|_ e a+b(N-Z) ’ (5) 


with values Vl = 160.0 MeV, Vd = 15.0 MeV, a = 7.2 
and b = —0.3 adjusted to obtain the best possible de¬ 
scription of available half-life data [36]. 



FIG. 1. (color online) /3-decay Q-values for the silver and cad¬ 
mium isotopic chains. Values calculated in this work, marked 
with red squares, are compared with the values obtained from 
the FRDM mass table [37], denoted with blue diamonds, the 
HFB-21 model [38] and with existing experimental data [36] 
denoted with black circles. 

The model employed in this study is well suited to 
describe the properties of the ground state and of the 
excited states of even-even nuclei. However, as the r- 
process naturally includes odd-A and odd-odd nuclei as 
well, it becomes necessary to provide at least an approx¬ 
imate description of the decay properties of these nuclei. 
With this aim, we compute the ground state of odd nuclei 
by employing the same model as for even-even nuclei, but 
constraining the expectation value of the particle number 
operator to an odd number of protons and/or neutrons. 
In this way an even RHB state is obtained, whose energy 
is different from true odd nucleus ground state energy by 
energy of the odd quasi-particle [39, 40]. To quantify 
the validity of this approximation we show the computed 
Qp values for the silver and cadmium (Z = 45) isotopic 
chains in Fig. 1, calculated by approximating the binding 
energy of the daughter nucleus [22]: 


Q — AM nH + A n — A p — -E 2? p 

= M ex (Z, N) — M ex (Z +1,N — 1), (6) 

where A„ and A p are the neutron and proton Fermi ener¬ 
gies in the ground state of the parent nucleus, E 2 q P are 
the lowest quasiparticle energies. A M u h = 0.782 MeV 
is the mass difference between a neutron and the hydro¬ 
gen atom, and M ex ( Z , N ) is the mass excess. In silver 
isotopes (odd Z), the Qp values are well reproduced for 
nuclei with an even number of neutrons, and only slightly 
underestimated for odd-odd nuclei. The predicted values 
closely follow the jump above the closed neutron shell, 
and begin to deviate from the FRDM values after the 
shell closure, while the HFB-21 results are found be¬ 
tween the other two models. For the cadmium isotopes 
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the model overestimates the Qp values in nuclei with an 
even number of neutrons in lighter isotopes, but above 
N = 78 the results reproduce the experimental values ac¬ 
curately. Again, the model very successfully reproduces 
the changes around the closed neutron shell. For heavier 
isotopes there are no available data, but all the models 
coincide in their predictions of the Qp values. We can 
conclude from these results that this approach to esti¬ 
mating the Q-values provides a good agreement with the 
available data and does not hinder our study of weak- 
interaction processes. 


B. Beta-decay half-lives 


In the present study of beta-decay half-lives we include 
both allowed (L = 0) and first-forbidden (L = 0,1) tran¬ 
sitions. The beta-decay rate for a transition between an 
initial and final nuclear state is equal to [41] 

In 2 f Po 

A = JfJ Q vl (W 0 - W) 2 F(Z , W)C(W)dp e , (7) 

where W is the electron energy in units of m e c 2 with 
Wo being the maximum electron energy that is equal 
to the difference in nuclear masses between initial and 
final nuclear states, Wo = (M* — Mf)/m e , and p e is the 
electron momentum in units of m e c. We approximate the 
maximum electron energy with [ 22 ]: 

Mi — Mf ~ A n — A p + A m np — Eqrpa • ( 8 ) 

The constant K is measured in superallowed beta-decay 
to be, K = 6144 ±2 s [42], 

The shape factor C(W) differs for various decays. For 
allowed decays it is energy independent. In the case of / 3~ 
decay of neutron-rich nuclei the shape factor is simply the 
Gamow-Teller reduced transition probability, C(W) = 
B(GT), with 


B(GT) = g\ 


</||Efc cr ^- 

(2 Ji + 1 ) 



(9) 


The matrix element is reduced with respect to the spin 
operator cr only using the Condon-Shortley phase convec¬ 
tion [43]. The sum runs over all nucleons. For the isospin 
lowering operator we use the convention t_ |n) = | p). Fi¬ 
nally, gA = —1.2701(25) [44] is the weak axial coupling 
constant. 

For first-forbidden transitions the shape factor reads: 


where the factors k, ka , kb and kc are defined as [45]: 

1 


k = 


+ 9 W 


- 1 ( 0 ) 


Cl +^( x + u) 2 - ^171 u(x + u) 


+ jg W o( 2x + u ) 2 - Jg X z( 2x ~ u ) 2 


(i) 


(lla) 


~ 2 (W 2 - A 2 ) 


( 2 ) 


ka = 


~\uY - ^Wo (4a ; 2 + 5m 2 ) 
o y 


-i(i) 


W 0 z 2 


- 1 ( 2 ) 


kb = jn i7i {- M( 0 ) + [Ci(z + «)](i)} , 
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kc = — [ 8 m 2 + (2a; + m ) 2 + A 2 (2a; — m) 2 ] 
18 


(i) 


cub) 

(lie) 

(lid) 


+ Y 2 K 1 + X ^ z2 ^(2) > 


( 2 ) ' 


with 


V = £'v + £w', Co = V -)- gMlWo, 

Y = €'y- C(u' + x'), Ci = Y + A(m - x)W 0 . 


( 12 ) 


The numbers in parenthesis after the closing brackets 
denote the rank of the operators inside. The parameter 
71 is given by \J\ — (aZ) 2 . For the coulomb wave func¬ 
tions we use the approximations ~ 1 and A 2 ps 1 [46]. 
The quantity £ = aZ/(2R) where we choose R = a/ (r ) 2 . 

In the Condon-Shortley phase convention [43] the ma¬ 
trix elements for j3~ transitions are: 


= -9aV?> 


(f\\E k r k [Ci <2> cr k ] ° t k 


x = — 


V2 Ji + 1 

(/ ||Efc ffcCi £-|| i) 


u = —g A Y2 


y/2Ji + 1 

(f\\Z k r k [Cfa**] 1 ** 


2 = 2 g A 


(/||E fc r fc [Cf 


V2 Ji + 1 


<cr k Vt k _ 


y/2Ji + 1 


(13a) 

(13b) 

(13c) 

(13d) 


w = -g A 


(f Efc r kI(Y 1,1,1, r k ) [C'l <g> a k ] ° t k _ i) 

X -W 
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u = -g A 


a /3 a /2 Ji + 1 

2(/||E fc r- fc /(l,l,l,l,r fe )C^||) 


a/ 2 Jj + 1 


(13f) 


2 a /2 (/ Efc r k I( 1 , 1 , 1 , 1 , r k ) [C k ® cr k ] 1 t k _ i) 

X ---(±4) 


The matrix elements connected with the relativistic cor¬ 
rections are 


Yv = -g A 

?v = - 


\/2Ji + 1 
with the ] 

</||E fc 75 fc £||») 

a/2 Ji + 1 

(/ ||a fc t(i || i) 


(14a) 


a/ 2 Ji + 1 ’ 


C(W) = k + kaW+ kb/W+ kcW 2 . (10) 


(14b) 
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with 75 and a being the Dirac matrices. Due to the fact 
that our formalism is relativistic we do not perform a 
non-relativistic reduction in the evaluation of the matrix 
elements in equations (14). The quantity Clm equals 

Clm = \/ 2TTi Y lM ' (15) 

where Ylm are the spherical harmonics 

As the matrix elements are independent of the elec¬ 
tron and neutrino energies, the integrals over the elec¬ 
tron phase space are evaluated independently and appear 
only as multiplicative factors. The function 1(1, 1,1,1, r) 
takes into account the nuclear charge distribution, and in 
the approximation of the uniform spherical distribution 
it has the form [45] 


/(l,l,l,l,r) 


3 

2 


1 -Hs ) 2 0 <r<i? 

f-Mf) r>R 


(16) 


Systematic calculations of beta-decay half-lives deter¬ 
mined by Gamow-Teller transitions have shown that the 
theoretical matrix elements need to be quenched by a 
factor q that has been found to be independent of the 
particular transition considered and approximately con¬ 
stant over the nuclear chart [47, 48]. There is also ev¬ 
idence that the total first-forbidden transition strength 
is overestimated by theoretical approaches. In addi¬ 
tion, the quenching factor seems to depend on the model 
space used and on the particular first-forbidden opera¬ 
tor [18, 49-53]. So far no study has addressed the quench¬ 
ing of first-forbidden transitions of global calculations of 
beta-decay half-lives. In our calculations, we account 
for the quenching of both the Gamow-Teller and first- 
forbidden transitions by using the same effective value of 
g A = -1.0. 


III. RESULTS 



neutron number neutron number 

FIG. 2. (color online) Comparison of the calculated half- 
lives with the FRDM results and the experimental data for 
the strontium (upper left), yttrium (lower left), silver (upper 
right) and cadmium (lower right) isotopic chains. In cases 
where they are not visible, the error bars are smaller than the 
data marker. 



FIG. 3. (color online) Percentage of first-forbidden transitions 
contributing to the decay of Sr, Y, Ag and Cd isotopes. In 
the cases where it was not possible to extract the value for a 
particular model, the data point was not plotted. 


A. /3-decay half-lives 

In Fig. 2 we compare the calculated half-lives for four 
isotopic chains, including both the even and odd atomic 
numbers, where the half-lives are obtained from the decay 
rate in Eq. (7): 


In the four panels we show a comparison of the present 
results with the data and the FRDM results for the stron¬ 
tium (Sr, Z = 38), yttrium (Y, Z = 39), silver (Ag, 
Z = 47), and cadmium (Cd, Z = 48) isotopic chains. 
The two lightest isotopic chains, strontium and yttrium, 
are found between two closed proton shells Z = 28 and 
Z = 50. Our calculation has difficulties reproducing the 
half-lives of the longest lived nuclei due to the rather low 


Q-values, however for isotopes more neutron-rich than 
N = 58 we obtain very good agreement with the data, 
with a slight tendency to underestimate the half-lives of 
yttrium isotopes. The FRDM+QRPA approach displays 
erratic jumps in the half-lives and odd-even staggering 
that is not present in the data. For very neutron-rich 
isotopes the two models provide very similar results due 
to the increase of Q-beta values. 

The right panels compares the present calculation and 
the the FRDM+QRPA model with data for nuclei close 
to the Z = 50 closed proton shell, including the silver and 
cadmium isotopic chains. Similar to the two lighter iso¬ 
topic chains, the D3C* interaction provides an excellent 
description of the decay properties of neutron-rich nu¬ 
clei, with the results correctly following the trend along 
the chain. In the case of silver, the present work re¬ 
produces the trend quite well in comparison with the 
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FRDM. Most likely, because the FRDM reproduces the 
half-life of 114 Ag it cannot provide a good trend for most 
of the isotopic chain. The half-life of this nucleus is un¬ 
characteristically short compared with the surrounding 
isotopes due to the J 71 ' = 1 + angular momentum and 
parity of the ground state which quickly decays into the 
0+ ground state of 114 Cd. For cadmium nuclei both mod¬ 
els do quite well, although the FRDM shows unexpected 
changes in the half-lives along the chain, as in the cases 
of other elements. 

In the following figure, Fig. 3, we plot a comparison 
of the contribution of the first-forbidden transitions in 
the total decay rate, in percent, between the two models. 
For the FRDM, we have extracted the values by com¬ 
paring the results of Ref. [54] and Ref. [19]. Due to the 
changes in the details of the calculation between the two 
publications, the values presented here are not precise. 
We expect them to be a good qualitative measure of the 
actual values. The FRDM calculation predicts, in two 
out of four isotopic chains, a large contribution (more 
than 50%) of first-forbidden transitions in nuclei close to 
stability, and a decrease in more neutron-rich nuclei. For 
nuclei in the isotopic chains of yttrium and silver, FRDM 
predicts a steady, small contribution below 20 % for most 
of the nuclei, with odd-even staggering. 

The present results predict a rather small contribution 
of the forbidden transitions in nuclei close to the stability, 
with a smooth increase with additional neutrons. This 
increase is particularly visible at the N = 82 neutron 
shell closures for Ag and Cd isotopes, where the neutrons 
begin to occupy the 2 / 7/2 and the I/ 19/2 orbits, enabling 
the strong u2 / 7 / 2 —> 7rl 57/2 and v\h§i 2 —> 7rlgg/ 2 transi¬ 
tions, which significantly increase the contribution of the 
first-forbidden transitions as these transitions appear at 
low excitation energies in the daughter nucleus. In the 
two lighter isotopic chains, the contribution is negligi¬ 
ble for the lightest isotopes and increases smoothly up 
to ~ 30% where it reaches a plateau. The main differ¬ 
ence in the transitions between e.g. 94 Sr and 114 Sr is the 
occupied neutron 2d^/ 2 state, which contributes to the 
total rate by decaying to the proton 2 p states. Above 
the N = 82 closed neutron shell, the contribution of the 
first-forbidden transitions again rises with the occupation 
of the 2 / 7/2 and the 1 /i 9 / 2 states, just as in the silver and 
cadmium isotopes. 

As discussed above, our calculations are only strictly 
applicable to even-even nuclei. However, we have ex¬ 
tended them to odd-A and odd-odd nuclei in order to 
provide a complete coverage for r-process simulations. 
Our approach cannot account for mismatches between 
the angular momenta of the ground states of parent and 
daughter nucleus in the decay. However, this aspect be¬ 
comes less and less important as the Q- value for the decay 
increases. This is one of the reasons why our approach 
becomes better as we move far from stability. Our predic¬ 
tions for half-lives are certainly not worse that those com¬ 
puted by the FRDM+QRPA approach [19] that is cur¬ 
rently the standard model for half-lives used in r-process 


simulations. We will discuss this aspect in more detail in 
the next section. 

It is also possible to compare the present results with 
the newest data available to assess the quality of the 
extrapolation of half-lives into currently inaccessible re¬ 
gions. Fig. 4 compares our calculations with the re- 



FIG. 4. (color online) The ratio of the calculated and mea¬ 
sured half-lives for nuclei in the 36 < Z < 43 isotopic chains 
from Ref. [9], for D3C* interaction (full line) and the FRDM 
(dashed line). Different colors denote the isotopic chains, 
while the thick black lines denote deviations by a factor 2. 

cently measured beta-decay half-lives [9] for a number of 
isotopic chains with atomic number in the range 36 < 
Z < 43. The results obtained with the D3C* interaction 
are in agreement with the experiment within a factor 2 . 
The FRDM+QRPA approach, on the other hand, tends 
to predict too long half-lives and displays jumps in the 
ratio for specific nuclei. 



FIG. 5. (color online) Comparison of the data (denoted with 
black circles) with the results of the present work (red squares) 
and the FRDM (blue diamonds). 

As we move to heavier nuclei both models have signifi¬ 
cant difficulties while attempting to reproduce the recent 
data. Fig. 5 compares the two models with measurements 
[12, 13]. FRDM significantly overestimates the half-lives 
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of most of the nuclei, by more than an order of mag¬ 
nitude. It is only for nuclei around the Z = 82 closed 
proton shell that provides a reasonable description of the 
data. An interesting effect occurs in the case of rhenium 
isotopes, where the half-life of 194 Re is overestimated, 
but half-lives of the two heavier isotopes are reproduced. 
The discrepancy seems to be related to the change of pro¬ 
late to oblate deformation in moving from 194 Re to the 
heavier isotopes [37]. With the exception of 194 Rc the 
D3C* calculations reproduce all beta-decay half-lives of 
the Z < 79 nuclei shown in figure 5 within an order of 
magnitude. Larger differences appear for heavier nuclei. 
Taking into account that the nuclei shown in figure 5 have 
half-lives in the range 1 s to 100 s one should conclude 
that both FRDM and D3C* approaches provide a fair 
reproduction of data. Certainly, the agreement is worse 
that the one obtained in figure 4 but this should not be 
consider as a drawback of the two models discussed to re¬ 
produce half-lives of heavier nuclei. It simply reflects the 
fact that for heavy nuclei the half-lives are only known 
for relatively long lived nuclei close to the stability that 
normally have small Q-beta decay windows. The theo¬ 
retical description of these decays is a challenge for any 
theoretical model as the decay rate depends on very few 
selected transitions. 
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with the very recent measurement of half-lives of 110 
neutron-rich isotopes, 40 of which have not been previ¬ 
ously measured. In the upper panel we plot the results of 
the present study, from the isotopes of rubidium to tin, 
and from N = 65 to N = 89. Almost all of the values 
can be found within a factor of 2 from the data, with 
very little scatter for different isotopic chains, except for 
the case of tin where we underestimate the half-lives by 
a factor of 5. Within a particular chain the data is re¬ 
produced smoothly, with very weak odd-even staggering. 
In the case of FRDM the scatter is much larger, and we 
also observe some regularities. The half-lives of nuclei in 
isotopic chains with an even atomic number are system¬ 
atically overestimated by a factor of two or more - these 
are the Zr, Mo, Ru, Pd, Cd and Sn isotopes. 
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FIG. 7. (color online) The logarithm of the ratio of the calcu¬ 
lated and experimental half-lives is plotted for all neutron-rich 
nuclei with available data, with half-lives shorter than 1 hour. 
Black squares indicate that the model did not provide a finite 
half-life. In the top panel we plot the results of the present 
study, while in the bottom panel we plot the ratios obtained 
with the FRDM. 


FIG. 6. (color online) The ratio of the calculated and experi¬ 
mental half-lives for the recently measured 110 isotopes [55]. 
In the top panel we present the results obtained in the present 
study, while in the bottom panel we present the FRDM re¬ 
sults [19] . 

Finally, in Fig. 6 we compare the theoretical results 


Fig. 7 shows the ratio of the calculated and measured 
half-lives for neutron-rich nuclei with half-lives shorter 
than 1 hour. We show the results for the FRDM+QRPA 
approach and our own calculations (bottom and top pan¬ 
els, respectively). The present calculation has difficulties 
in reproducing the half-lives for nuclei close to simultane- 















































ous neutron and proton shell closure. In particular, the 
greatest discrepancies appear in the regions “north-west” 
of the doubly closed nuclei 78 Ni and 132 Sn. The amount 
of low energy GT strength for such nuclei is very sensitive 
to correlations [18] that are not captured by the QRPA 
approach. For an improved description of the low-energy 
strength distributions of these nuclei, complex configura¬ 
tions need to be included to shift the transition strength 
from the GT resonance to low energies [56, 57]. Cou¬ 
pling to collective phonon excitations is known to have 
a strong impact on the single particle structure around 
shell closures [58], and will enhance the decay rate of 
these nuclei [59]. 

There are also three regions where the D3C* calcula¬ 
tions predicts shorter decay half-lives, i.e. the regions 
around the 90 Kr (Z = 36, N = 54) and 140 Ba (Z = 56, 
N = 84) nuclei, and the Z <82 N > 126 region. 
The first two cases correspond to regions of low defor¬ 
mation [60]. However, there are indications that nuclei 
on that region are gamma soft [61, 62]. So far the im¬ 
pact of gamma-softness on beta-decay rates has not been 
explored. It is interesting that the FRDM+QRPA ap¬ 
proach also predicts shorter half-lives for these two re¬ 
gions. In the case of the region Z < 82 N > 126, we 
are dealing with nuclei with half-lives longer than 100 s. 
The decay rate is determined by few transitions that 
depend very sensitively on the predicted single-particle 
structure. In particular, for Hg and T1 isotopes, above 
the N = 126 closed neutron shell the decay is domi¬ 
nated by the Gamow-Teller vlin /2 —> ^ 1 * 13/2 an( i the 
1“ 1 ^ 1 * 11/2 —> Trlhg /2 configurations. These transitions 
are not possible for isotopes with N < 126, but above 
the neutron shell closure the vi\\i 2 is becoming signifi¬ 
cantly occupied (occupation probability is already ~ 0.2 
for 210 Hg and 211 T1) at an energy of E ill/2 « —2.5 MeV. 

Fig. 8 shows the ratio of the calculated and the exper¬ 
imental half-lives versus the experimental half-life. Dif¬ 
ferent symbols are used for even-even, the odd-N, odd-Z 
and the odd-odd nuclei. The data is taken from Ref. [36]. 
For long lived nuclei we observe a very large spread of the 
ratios due to the high sensitivity of the decay rate to the 
details of the low energy transition strength. Short-lived 
nuclei have large Q-values and generally several transi¬ 
tions contribute to the decay rate. This makes these nu¬ 
clei less sensitive to the detailed position and strength of 
particular transitions. Our results show that the ratio of 
half-lives converges towards one with decreasing half-live, 
for even-even, odd-A and odd-odd nuclei equally. This 
behaviour points to the reliability of the model for all 
types of nuclei taking part in the r-process, as almost all 
of the predictions lie within an order of magnitude from 
the data, with a large majority within a factor of 2. 

In Table I we provide a quantitative analysis of the 
results, by using several measures of the model precision. 
Because both the observable (half-life) and the ratio of 
the calculated and experimental values span orders of 
magnitude we define rt as the logarithm of the ratio of 
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FIG. 8. (color online) The ratio of the calculated and the 
experimental half-lives [36] versus the experimental half-life. 
Different types of nuclei are denoted with different symbols. 
Even-even nuclei are denoted with blue squares, odd-N with 
orange triangles, odd-Z with black diamonds, and odd-odd 
nuclei are denoted with red circles. 


calculated and experimental half-lives 

^calc. 

^=log 10 7^7, (18) 

4 1/2 

with the average value of the ratio providing a measure 
of the global deviation from experimental values 

1 N 

( 19 ) 

i=l 

while the standard deviation provides a measure of the 
spread [63] 



i=l 


The average deviation and the spread obtained using the 
D3C* and the FRDM model [19] are shown. The top 
part of the table provides data for nuclei with the ex¬ 
perimental half-life shorter than a particular value. The 
average value of deviation from experiment r is compara¬ 
ble between the two models. There are small differences 
for particular sets of nuclei, but they are not significant. 
The resulting standard deviation is larger with the D3C* 
for long-lived nuclei, confirming that the description of 
nuclei with small Qp values is challenging for the model. 
For nuclei with half-lives shorter than 1 s, the models are 
comparable both in the average and the standard devia¬ 
tion, and for even less stable nuclei D3C* even manages 
to provide a smaller spread of ratios. 










9 


TABLE I. Comparison of the average deviations (Eq. (19)) and their standard deviations (Eq. (20)) between the model 
used in the present work and the FRDM. In the upper part of the table, we show the results structured with respect to the 
experimental half-life of the nuclei included, while in the lower part of the table we show the results for nuclei with half-lives 
shorter than 1 s, but separated into groups of even-even, odd-Z, odd-N and odd-odd nuclei. 



r 

D3C* 

a 

r 

FRDM 

a 

Tex p. 


< 1000 s 

0.011 


0.889 

0.021 


0.660 

< 100 s 

0.057 


0.791 

0.040 


0.580 

< 10 s 

0.061 


0.645 

0.046 


0.515 

< 1 s 

0.011 


0.436 

0.019 


0.409 

< 0.1 s 

0.041 


0.195 

0.021 


0.354 

nucleus type 


even-even 

-0.037 


0.331 

0.333 


0.226 

odd-Z 

0.054 


0.328 

-0.128 


0.288 

odd-N 

-0.086 


0.387 

0.124 


0.436 

odd-odd 

0.089 


0.582 

-0.179 


0.409 

total 

0.011 


0.436 

0.019 


0.409 


In the bottom part of the table we provide the results, 
including nuclei with the experimental half-life lower or 
equal to 1 s, but shown for each type of nuclei, i.e. even- 
even, odd-Z, odd-N and odd-odd, separately. From this 
point of view, the results of the two models are sig¬ 
nificantly different. In the present work, the half-lives 
of even-even and odd-N nuclei are somewhat underesti¬ 
mated by roughly 10% - 20%, and the half-lives of odd-Z 
and odd-odd nuclei are, on average, overestimated by 
15% and 20%, respectively. These results agree with the 
assumption that the model is capable of describing all 
types of nuclei equally well, as the half-lives of odd-A 
and odd-odd nuclei are reproduced as well as half-lives 
of even-even nuclei. The difference between various types 
of nuclei becomes more evident by looking at the stan¬ 
dard deviations. Here, the spread of values for even-even 
and odd-A nuclei is comparable, but the standard devi¬ 
ation of the results for odd-odd nuclei is almost twice as 
large, pointing to potential angular momentum issues in 
the ground state of an odd-odd nucleus. It is interesting 
that our model does particularly well for even-N nuclei 
that are the most relevant nuclei for r-process nucleosyn¬ 
thesis as these nuclei are favored by neutron-captures. 

The FRDM+QRPA half-lives show a distinctly dif¬ 
ferent behaviour. As opposed to the D3C* model, the 
FRDM overestimates the half-lives of even-even and odd- 
N nuclei, by a factor of 2.15 and 1.33, respectively. That 
the predicted half-lives of even-even nuclei are on average 
twice the experimental value is surprising as, similarly to 
our calculation, the model should be more suitable for 
these nuclei. This fact may have important consequences 
for r-process nucleosynthesis and is might be responsible 
for the differences in r-process nucleosynthesis discussed 
on section HID. However, waiting point nuclei also in¬ 
clude odd-Z, even-N nuclei, and in both models there is 
a cancellation between even-even and odd-Z nuclei. 

On the other hand, the half-lives of odd-Z and odd- 
odd nuclei are underestimated by approximately 30% and 


50%, respectively. The standard deviations are compa¬ 
rable between the two models for odd-A nuclei, but are 
smaller for even-even and especially for odd-odd nuclei. 
Still, for both models the odd-N and odd-odd nuclei are 
the most difficult to describe and result with the largest 
spreads. So, even though the average f of the FRDM 
model is close to zero, it is actually a result of the can¬ 
cellation of overestimated and underestimated half-lives 
in different types of nuclei. This behaviour has also ap¬ 
peared in previous evaluations of the /3-decay rates (see 
Ref. [54], Table B.), where the half-lives of short-lived 
even-even nuclei were, on average, overestimated by al¬ 
most a factor of 4, while the half-lives of odd-A and odd- 
odd nuclei were underestimated by a factor of 2. In com¬ 
parison, the results of the present work are more con¬ 
sistently close to reproducing the data, even though the 
larger standard deviations highlight some unresolved is¬ 
sues. 



FIG. 9. (color online) The ratio of halflives obtained by the 
two models. 

Finally, in Fig. 9 we plot the ratio of half-lives obtained 
in the present study with the FRDM half-lives [19]. For 
light and medium-heavy nuclei, the results are quite simi¬ 
lar. Close to the valley of stability, the D3C* model tends 
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to provide longer half-lives than the FRDM for particu¬ 
lar nuclei, especially so in the regions “north-west” of 
the doubly closed nuclei 78 Ni and 132 Sn. Further from 
the valley of stability the two models provide compara¬ 
ble results, with the difference that the present calcula¬ 
tion predicts smoother increase in the decay rates. For 
very heavy nuclei with N > 126 however, the present 
study predicts half-lives to be shorted by more than an 
order of magnitude than the FRDM predictions. This is 
especially clear in nuclei with high atomic numbers, i.e. 
Z 95 where the difference is as large as three orders of 
magnitude. This may have significant consequences on 
the dynamics of the r-process nucleosynthesis in neutron 
star merger conditions [64] 


B. Impact of First-forbidden transitions 

The results presented in this manuscript are the first 
global calculations of beta-decay half-lives for neutron- 
rich nuclei that treat Gamow-Teller and first-forbidden 
transitions on an equal footing. In this section, we pro¬ 
vide a detailed analysis on the impact of first-forbidden 
transitions with particular focus on nuclei with magic 
neutron numbers N = 50, 82 and 126 that have been the 
subject of many theoretical calculations. Additionally, 
for the N = 50 and 82 experimental data is available 
that has contributed to constrain the calculations. 

Fig. 10 compares the beta-decay half-lives for r-process 
nuclei with neutron magic number N = 50,82 and 126 
with the FRDM+QRPA model [19], the shell-model cal¬ 
culations of ref. [18] and data [36]. For the N = 50 iso¬ 
tones (upper panel) with Z > 28 the present calculations 
fail to reproduce the measured half-lives. In these nuclei 
the f'l/ 2 proton shell is fully occupied suppressing low en¬ 
ergy Gamow-Teller transitions. The decay consequently 
proceeds mainly by first-forbidden transitions resulting 
in long decay half-lives. As a consequence the model 
predicts a large contribution of forbidden transitions to 
the decay rate (see fig. 11). This is probably a limitation 
of the QRPA type calculations that are not able to pro¬ 
duce enough correlations around the proton magic num¬ 
ber Z = 28. This problem does not seem to be present 
in the FRDM+QRPA approach. However, one should 
kept in mind that in this approach the Gamow-Teller 
and first-forbidden contributions are not derived based 
on the same microscopic model. As discussed above the 
FRDM+QRPA approach tends to produce strong odd- 
even effects that are not present in the data nor in the 
other approaches. 

The situation is different for the N = 82 and N = 126 
isotonic chains. In these nuclei there is no proton shell 
closure that hinders low energy Gamow-Teller transi¬ 
tions and in principle both Gamow-Teller and forbidden 
transitions are possible. The present study predicts a 
very similar contribution of first-forbidden transitions to 
the shell-model, ~ 20%, for all the N = 82 isotones. 
FRDM+QRPA predicts rather large forbidden contribu- 
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proton number 

FIG. 10. (color online) Beta-decay half-lives for N = 50 
(a), N = 82 (b), and N = 126 (c) isotopic chains. The 
present results are compared with those of FRDM [19], the 
shell-model [18], DF3 [27] and data [36]. 


tions for Z = 46,47,48 and similar values for the other 
isotones. The odd-even staggering present in the beta- 
decay half-lives is also manifested on the contribution 
of forbidden transitions. For N = 126 all models pre¬ 
dict an enhanced contribution of forbidden transitions. 
Ref. [18] argues that the contribution of forbidden tran¬ 
sition should decrease with decreasing proton number. 
This reduces the number of protons in the h\\i 2 increas¬ 
ing the role of ^/ig /2 —> 7r hn/ 2 Gamow-Teller transitions. 
While all models predict such a reduction, the magnitude 
varies from model to model. The present calculations 
predict a relatively minor reduction from 90% at Z = 73 
to 70% at Z = 66, while the shell-model predicts a much 
larger reduction reaching below 20% at Z = 66. Exper¬ 
imental data will be very important to further constrain 
the theoretical calculations. 

Fig 12 shows the contribution of beta-decay half-lives 
for all nuclei with predicted beta-decay half-lives shorter 
than 1 s. In general, first-forbidden transitions do not 
noticeably impact the total decay rate in nuclei close to 
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FIG. 11. (color online) The contribution of the first-forbidden 
transitions in the total decay rate for the N = 50, N = 82 and 
N = 126 isotonic chains, in the top, central and bottom pan¬ 
els, respectively. Present results, denoted with black squares, 
are compared with the FRDM calculation [19] (denoted with 
green triangles), DF3 calculation [20] (denoted with blue di¬ 
amonds) and a recent shell model calculation [18] (denoted 
with red circles). 
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FIG. 12. (color online) The contribution of the first-forbidden 
transitions in the total decay rate for nuclei with predicted 
half-lives shorter than 1 s. Black lines indicate the position 
of (predicted) closed proton and neutron shells. 


ity and consequently Gamow-Teller and forbidden tran¬ 
sitions contribute in similar amounts to the decay. The 
situation is different in the region between N = 50 and 
N = 82. Here, once protons completely fill the pf shell 
(Z > 40) forbidden transitions are suppressed and the 
decay is dominated by Gamow-Teller transitions. 

An interesting phenomena occurs once we move to nu¬ 
clei with N > 126 with a clearer manifestation once we 
reach N = 184. For these nuclei, neutrons and protons 
occupy orbits that differ in two units of the main oscil¬ 
lator quantum number, M = 2 n + l. Under the assump¬ 
tion of isospin symmetry the single particle orbits for 
neutrons and protons will be identical and consequently 
Gamow-Teller transitions will be exactly zero. However, 
due to isospin breaking mainly due to the coulomb inter¬ 
action the proton and neutron single particle states are 
not identical allowing for Gamow-Teller transitions, even 
if somewhat suppressed. This is the reason why forbid¬ 
den transitions dominate in this whole region with GT 
transitions representing around 20% of the decay rate. 


C. /^-delayed neutron emission 


the valley of stability with Z < 50. Here, protons and 
neutrons fill orbits with the same parity and the decay is 
dominated by Gamow-Teller transitions. It is interesting 
to notice the similar behaviour for nuclei located around 
' 8 Ni and 132 Sn. Moving in the “south-west” direction 
(Z < 28, N < 50 and Z < 50, N < 82 respectively) 
forbidden transitions are initially suppressed. However, 
once protons start to occupy orbits with different parity 
(Z < 20 or Z < 40) the contribution of forbidden tran¬ 
sitions grows to values of 40% to 60% depending on the 
detailed nuclear structure. Moving in the “east” direction 
(N > 50 and N > 82 respectively), we see that the con¬ 
tribution of forbidden transitions is rather independent 
of the charge number and does not vary when crossing 
the proton shell closure at Z = 28 and Z = 50, respec¬ 
tively. In these regions, for neutrons and protons there 
are valence orbits with both positive and negative par- 


/3-delayed neutron emission is another important com¬ 
ponent in the late stages of r-process nucleosynthesis. 
Here we approximate the probability of emission of x neu¬ 
trons as the ratio of the rates between S xn and S( x +u n 
separation energies to the total decay rate, i.e. 

min(Q /3 ,5( a , + i)„) 

i,Ei=Sxn \ 


Neutron emission probabilities are, therefore, very sensi¬ 
tive to the neutron separation energies and display strong 
odd-even effects. As mean-field models cannot repro¬ 
duce this staggering we employed the global nuclear mass 
model [60] to obtain the neutron separation energies. Ad¬ 
ditionally, only for the purpose of calculating the emis- 
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sion probabilities, we smeared the beta strength distribu¬ 
tions with a Lorentzian in order to reduce the sensitiv¬ 
ity of the neutron emission probabilities on the positions 
of particular transitions. The width of the Lorentzian 
was determined by calculating the average number of de¬ 
layed neutrons emitted by the fission fragments of 325 U, 
(n) = 0.0158 ±0.0005 [65]. In this way the width was ad¬ 
justed to T = 65 keV. From the probabilities, we obtain 
the average number of emitted neutrons after the decay 
of a nucleus as 


(n) = Y, iP in- ( 22 ) 

i 


In Fig. 13 we plot the average number of emitted neu¬ 
trons for neutron-rich nuclei included in this study. Close 
to the valley of stability the Q-values are smaller than 
the one neutron separation energies, making it impos¬ 
sible for the daughter nucleus to de-excite via neutron 
emission. With additional neutrons both the Q-value in¬ 
creases and the separation energies decrease enabling the 
emission of one, two or more neutrons after decay. This 
is particularly evident above neutron shell closures where 
the average number of emitted neutrons increases due to 
a drop in the separation energies. Another characteristic 
feature of (n) is the odd-even staggering as a function 
of atomic number. Isotopic chains with an even num¬ 
ber of protons typically have a smaller average number 
of emitted neutrons, i.e. smaller probabilities of emit¬ 
ting more neutrons, than their odd-proton neighbouring 
chains which mimics the behaviour of the Q-values. A 
smaller odd-even staggering is also present within any 
isotopic chain as a consequence of the staggering of the 
neutron separation energies. 

While both the D3C* and the FRDM models repro¬ 
duce these general features, in other aspects there are 
significant differences. Firstly, the FRDM evaluation 
was limited to a maximum of 3 emitted neutrons per 
decay, while the present study takes into account the 
possibility for up to 5 neutrons to be emitted in a de¬ 
cay. The consequence of this change is evident in very 
neutron-rich nuclei where the neutron separation ener¬ 
gies are very small, and where the D3C* model predicts 
a larger average number of emitted neutrons which con¬ 
tinuously increases up to the drip-line, while the FRDM 
model predicts a saturation of (n). On the other hand, 
for isotopes closer to stability the D3C* model predicts 
smaller vales of (n) and a more gradual increase towards 
the neutron drip-line. Finally, in very neutron-rich nuclei 
the FRDM displays decidedly unphysical, sudden oscil¬ 
lations between the values of (n) = 1 and (n) = 2 while 
the present study predicts a smooth increase with the 
neutron number. 
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FIG. 13. (color online) Average number of emitted neutrons 
after /3-decay calculated with the FRDM and the D3C* mod¬ 
els. Note that because the available FRDM results only con¬ 
tain the emission probabilities for up to 3 neutrons emitted, 
and in the present work the probabilities are obtained for up 
to 5 emitted neutrons, the colour scales are different. 


D. Consequences for heavy-element 
nucleosynthesis 


To explore the impact of the obtained results on heavy 
element nucleosynthesis, wc have performed r-process nu¬ 
cleosynthesis calculations based on the “hot” and “cold” 
r-process conditions from ref. [3]. In the “hot” r-process 
scenario the r-process evolves under an (n, 7 ) «=^ ( 7 , n) 
equilibrium that breaks down once the r-process freeze- 
out is reached. For the “cold” r-process scenario the 
evolutions proceeds by a competition between neutron- 
captures and beta-decays. The neutron capture and pho¬ 
todissociation rates used on the network calculations are 
based on the statistical model approach [ 66 ] using the 
FRDM mass model [37]. The only difference between the 
calculations is the use of the FRDM±QRPA beta-decay 
half-lives and beta-delayed neutron emission probabili¬ 
ties [19] or those computed in the present study. 

In Fig. 14 we compare the abundances obtained with 
the new decay half-lives and delayed neutron emission 
probabilities with the results obtained using the FRDM 
values for both the hot and cold r-process trajectories. 
One can notice important differences particularly in the 
region of the third r-process peak, A ~ 195 between the 
calculations. There are also differences in the abundances 
around the second peak but the astrophysical conditions 
explored by the present trajectories are expected to con- 
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FIG. 14. (color online) The abundances of heavy nuclei re¬ 
sulting from a hot r-process (upper panel) or a cold r-process 
(lower panel) calculation using the half-lives obtained from 
the FRDM and the current model. The solid circles denote 
the solar r-process abundances. 


tribute mainly to the third r-process peak region, see 
ref. [3]. To get further insight in the origin of the differ¬ 
ences in the abundances we show in figure 15 the evolu¬ 
tion of the neutron-to-seed ratio in the present calcula¬ 
tions. We clearly see that for those calculations that use 
the D3C* the r-process freeze-out, defined as the moment 
when the neutron-to-seed ratio reaches a value of 1, oc¬ 
curs at an earlier time. An earlier freeze-out has impor¬ 
tant consequences for the shaping of the final abundance 
pattern during the decay to stability. This is because it 
occurs at a higher value of the density that speeds up 
the freeze-out producing a faster and larger drop in the 
neutron-to-seed ratio. As a consequence, there is less 
time for neutron-captures during the freeze-out. These 
neutron captures are in fact responsible for shifting the 
peaks to higher A values as discussed in ref. [3]. 

The evolution proceeds faster with the new half-lives 
because they tend to be shorter than those computed 
by the FRDM+QRPA approach. This is particularly 
the case near to the magic shell closures N = 82 and 
N = 126 where the r-process moves closer to the stabil¬ 
ity and reaches the longest living nuclei. Both for the hot 
and cold r-process conditions, the r-process has reached 



FIG. 15. (color online) The evolution of the neutron-to-seed 
ratio during the r-process. The upper panel shows the results 
for the hot r-process while the lower panel for the cold r- 
process. Dashed lines denote the results obtained with the 
FRDM decay data, and full lines denote results obtained in 
the present study. 


beta-flow equilibrium by the freeze-out time. Under these 
conditions the product of the average beta-decay rate 
times the abundance for each isotopic chain is constant, 
i.e. X /3 (Z)Y(Z) = constant, with Y(Z) = ^ 4 y(Z, A) 
and 

xp{z) = ^Z) = Y(Z) 52 X f>( Z > A)Y(Z, A). (23) 

This means that the isotopic abundances, Y(Z), are pro¬ 
portional to the lifetime, t(Z ), and the material accumu¬ 
lates in regions with the longest lifetimes that correspond 
to the magic shell closures. Figure 16, shows the isotopic 
lifetimes at the freeze-out for the different calculations. 
One can clearly see that the lifetimes are substantially 
longer for the FRDM+QRPA approach particularly close 
to the magic neutron numbers N = 82, corresponding to 
Z ~ 50, and N = 126, corresponding to Z ~ 70. 































14 



FIG. 16. (color online) The isotopic lifetimes defined in Eq. 
(23) as a function of the atomic number, at the moment of 
r-process freze-out. In the top panel we show the results for 
the hot r-process, and in the bottom panel the results for the 
cold r-process. 


We expect that the variations discussed above will re¬ 
main for a broad range of r-process conditions and in 
particular also for neutron-star merger conditions. This 
aspect will be explored in more detail in a future pub¬ 
lication. An important aspect for r-process calculations 
under r-process conditions is the energy generation dur¬ 
ing the r-process. This depends mainly on the beta-decay 
rates and determines the evolution of temperature dur¬ 
ing the r-process [67]. The nuclear energy production 
is also important to determine the evolution of material 
ejected in highly eccentric loosely bound orbits. Depend¬ 
ing on the energy produced by the r-process this mate¬ 
rial can become unbound affecting the late time evolu¬ 
tion of short gamma-ray burst afterglows assuming they 
are powered by matter fallback [68]. At later times the 
radioactive heating from the beta-decay of r-process ma¬ 
terial will also affect the neutron-star merger remnant 
evolution [69] and the electromagnetic transient, kilo- 
nova, lightcurve [70-72]. For all these applications it is 
important to know how much of the energy produced 
by beta-decay is lost in the form of neutrinos and how 
the remaining energy is distributed between the electrons 
and gamma-rays emitted from nuclear deexcitation. We 
have determined for each decay these quantities based on 
the computed beta-decay strengths. A small part of the 
energy will be carried away by the neutrons emitted af¬ 
ter the beta-decay. Our gamma channel does in fact also 
include this energy as it is assumed that the emitted neu¬ 
trons will thermalize in a relatively short timescale. All 
this information is included in the supplemental material 
that also includes the relevant beta-decay rates. 

Figure 17 shows the average energy of electrons, neu¬ 
trinos and photons produced after beta-decay as a func¬ 
tion of the decay Q- value. The red lines denote the one- 
half and one-third ratios between the energy and the Q- 





Q [MeV] 


FIG. 17. (color online) Average energy of the electrons, neu¬ 
trinos and photons produced after beta-decay 


values. 


4V, 7 (Q) = \q, El^iQ) = \q- ( 24 ) 

One can see that the ratio of the average electron en¬ 
ergy and the nuclear Q -value is concentrated around one 
third, with the same artio for the antineutrinos is in the 
range 0.3-0.5. r-process nucleosynthesis is not expected 
to be sensitive to this range of variations in the heating 
rate [67]. However, kilonova light curves are more sensi¬ 
tive to this range of variation [70] and our calculations are 
the first to provide this information for future modeling. 


IV. CONCLUSION AND OUTLOOK 

In this study we have performed a large scale calcu¬ 
lation of /3-decay rates and ^-delayed neutron emission 
probabilities for neutron-rich nuclei with 8 < Z < 110. A 
fully self-consistent theoretical framework was employed 
based on the relativistic nuclear energy density func¬ 
tional. The advantages of this approach are: (i) self- 
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consistent modeling of all relevant transition matrix el¬ 
ements without any adjustment of the model parame¬ 
ters to the nucleus under consideration, and (ii) equal 
treatment of allowed and first-forbidden transitions. The 
r-process involves even-even nuclei, but also odd-A and 
odd-nuclei. It is, therefore, very important that all types 
of nuclei are treated within a single theoretical frame¬ 
work. We have tested the quality of the description of 
odd nuclei by comparing the /3-decay Q -values with those 
of the FRDM and the data for the silver and cadmium 
isotopic chains, where the results show very good agree¬ 
ment with the available data. Additionally, in agreement 
with data there is no noticeable odd-even staggering in 
the obtained half-lives, and the statistical analysis shows 
only a small tendency to overestimate the half-lives of 
odd-Z and odd-odd nuclei. Additionally, even though 
the model used assumes a spherical shape of the nucleus 
under consideration, it has been shown that it can ac¬ 
curately describe the reaction rates both in the regions 
of the nuclear chart close to closed shells and in open- 
shell nuclei known to be well-deformed. Thus, the model 
employed is shown to reproduce /3-decay data across the 
whole neutron-rich side of the nuclear chart with good 
precision, regardless of the type of nucleus. 

In general, the key features of the decay rates across 
the nuclear chart have been reproduced. We obtain a 
very good agreement with the data for isotopic chains 
with both an even and odd atomic number, without 
any noticeable odd-even staggering along the chain and 
closely following the trend of the measurements. This 
characteristic ensures a smoother behaviour during the fi¬ 
nal decay towards the line of stability in the late stages of 
the heavy element nucleosynthesis. The FRDM results, 
however, often displays erratic jumps in the half-lives 
without an underlying physical reason, as is evident in 
e.g. strontium and cadmium isotopic chains. In compar¬ 
ison to recent data, both models produce similar discrep¬ 
ancies. In the light nuclei region, both models reproduce 
the data within a factor of 2, for a wide selection of nuclei 
with atomic numbers ranging from Z = 36 to Z = 43. In 
the region of very heavy nuclei, the FRDM overestimates 
the half-lives of nuclei below the Z = 82 closed proton 
shell, while the present work provides a reasonable de¬ 
scription of the decay properties. However, around the 
N = 82 shell, the D3C* model significantly underesti¬ 
mates the half-lives. Extending the analysis to the whole 
neutron-rich part of the nuclear chart, we obtain very 
good agreement with the data, except in the regions just 
above a closed proton shell, but below the next closed 
neutron shell, i.e. the same regions where the model pre¬ 
dicts high contribution of the first-forbidden transitions. 
We also performed a statistical analysis, where the stud¬ 
ied the quality of the half-life predictions in contrast to 
the measurements, and depending on the type of nuclei. 
We obtain a very good average reproduction of the data, 
with a tendency to overestimate the half-lives of odd-Z 
nuclei. In contrast with this result, the very good overall 
agreement of the FRDM predictions with the measure¬ 


ments is a consequence of a cancelation of a significant 
overestimate of the half-lives of even-Z (especially even- 
even) nuclei with a smaller underestimate of the half-lives 
of odd-Z nuclei. In both cases, the description improves 
for shorter lived nuclei, due to the importance of the 
lepton phase space that increases sharply with larger de¬ 
values. 

The consistent inclusion of the first-forbidden transi¬ 
tions in this study has allowed for a systematic study of 
the contribution of parity changing transitions to the de¬ 
cay rate for a large number of nuclei. In specific cases, 
namely the N = 50, N = 82 and N = 126 isotonic chains, 
the results of the present study agree well with date and 
results based on other models, except in particular cases 
such as the nuclei above the Z = 28 proton shell where 
we predict a strong contribution of the forbidden tran¬ 
sitions. For N = 126 isotones, the model predicts an 
increase of the contribution of first forbidden transitions 
to the total decay rate, in agreement with the results 
obtained by other models. For nuclei with N < 126, we 
predict a smooth increase of the contribution of the parity 
changing transitions from the valley of stability towards 
the neutron drip-line, and a saturation at approximately 
40%-60% of the total rate. Nuclei with N > 126 show 
a suppression of Gamow-Teller transitions that became 
only possible by isospin breaking terms mainly due to 
coulomb. 

Closely related to /3-decay is the process of /3-delayed 
neutron emission, which also affects the late stages of the 
r-process by contributing to the neutron flux and the re¬ 
distribution of mass. The characteristic quantity is the 
average number of emitted neutrons which is a measure 
of the decay rate that falls above the neutron separation 
energy. This quantity depends on both the Q-value and 
the one, two and more neutron separation energies and 
displays odd-even staggering both between the isotopic 
chains and along a particular chain. We have extracted 
the probabilities for emission of one or more neutrons and 
the average number of emitted neutrons from the previ¬ 
ously calculated decay data. The results show a physical 
saturation of the average number of emitted neutrons 
closer to the neutron drip-line which indicates that the 
decaying matter could provide a significant number of 
neutrons. In contrast, the FRDM predicts a similar but 
overall lower number of emitted neutrons per decay, with 
oscillating behaviour close to the neutron drip line. 

The present results have been applied in an r-process 
calculation in a hot and cold r-process scenarios corre¬ 
sponding to matter ejected under moderately high en¬ 
tropy. We find substantial changes in the r-process abun¬ 
dances that are mainly due to the fact that for r-process 
nuclei the present half-lives are systematically shorter 
than those computed by the FRDM+QRPA approach. 

This study has shown that the theoretical framework 
based on the relativistic nuclear energy density functional 
is a mature one, capable of providing a good description 
of sensitive physical quantities. However, the description 
of half-lives of a large range of nuclei remains challenging, 


16 


especially for nuclei close to the valley of stability. Fur¬ 
ther advances in this field will require effort along three 
possible directions: (i) improving the relativistic energy 
density functional will serve to provide more accurate Qp 
values and enhance the description of the energy gener¬ 
ation during the decay released in the form of electron, 
antineutrinos and gamma-rays; (ii) a consistent treat¬ 
ment of deformations will be necessary to enhance the 
description of the transitions in particular open-shell nu¬ 
clei; and (iii) inclusion of complex-configurations in the 
QRPA basis will enrich the transition spectrum, espe¬ 
cially so in the low-energy region, and provide a better 
description of beta-decays, but also give information on 


the detailed structure of the transition strength. 
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